Microscopic mechanism for transient population inversion and optical gain in graphene 



Torben Winzerjj Ermin Malic, and Andreas Knorr 

Institut ftir Theoretische Physik, Nichtlineare Optik unci Quantenelektronik, 
Technische Universitat Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 

A transient femtosecond population inversion in graphene was recently reported by Li et al, Phys. Rev. 
Lett. 108, 167401 (2012). Based on a microscopic theory we clarify the underlying microscopic mechanism: 
Transient gain and population inversion in graphene occurs due to a complex interplay of strong optical pumping 
and carrier cooling that fills states close to the Dirac point giving rise to a relaxation bottleneck. The subsequent 
femtosecond decay of the optical gain is mainly driven by Coulomb-induced Auger recombination. 



Light amplification and opticalgain are of central impor- 
tance in the optical sciences (H U . They are based on the 
quantum effect of induced photon emission, and constitute 
one of the key ingredients in optoelectronic devices and op- 
tical data communication [2]. Therefore, there is an ongoing 
quest for gain materials in different frequency ranges. Particu- 
lar questions involve the mechanical stability of materials and 
the possibility of electrical injection. In recent years much 
research is devoted to graphene, a two dimensional sheet of 
carbon atoms 13-01. Beside its distinct transport character- 
istics ¥M, its optical properties have been intensively investi- 
gated ||5t l6Ll8l1 with respect to transient broadband and 
saturable f I UN 211 absorption. 

Recently, T. Li and coworkers 111 311 demonstrated also the 
possibility of optical gain based on population inversion in 
optically pumped graphene within a differential reflection 
measurement. The experiment reveals that gain appears at 
pump fiuences of approximately 2 mJ/cm 2 on a femtosecond 
timescale for optical transitions smaller than the pump energy 
| L3J] . Graphene as a coherent source of THz radiation was first 
suggested by Ryzhii et al. lfl4l [l5ll and experimental indica- 
tions were recently found in an optical non-differential mea- 
surement Il6ll . Since graphene has no band gap and its optical 
properties are dominated by distinct many body effects, the 
observation of optical gain in graphene is crucial to our un- 
derstanding of the light-matter interaction in this material. In 
particular, it is of outermost importance for graphene-based 
applications to understand the conditions for the occurrence 
of optical gain. 

In this Letter, we perform microscopic calculations on the 
carrier dynamics in the same excitation regime, as recently 
reported in experiments 11311 . confirming the possibility of 
gain and population inversion in graphene. We develop a 
complementary understanding with respect to the experiment, 
since our theory directly accesses the time- and momentum- 
resolved non-equilibrium carrier distribution during and after 
a strong optical excitation pulse. Therewith, we are able to 
analyze the predominant interplay of many body mechanisms 
(electron-phonon and electron-electron interaction) and light- 
matter interaction giving rise to the gain regime. In partic- 
ular, as the most important result of this paper, we identify 
as the crucial process the phonon-induced intraband scatter- 
ing, which fills optically active states at low energies. This 
process is responsible for a carrier relaxation bottleneck close 
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FIG. 1. (a) Scheme of optical pumping (red) and the build-up of a 
phonon bottleneck (black). Energetically lower probe pulse (blue) 
can be amplified via induced emission (yellow), if population inver- 
sion occurs, (b) Pump-induced change of occupation Apk (t) in the 
most efficiently excited state ko for different pump fiuences. We ob- 
tain a qualitatively and quantitatively excellent agreement with the 
experiment, cf. Fig. 2(b) in Ref. Hal , (c) In contrast to the exper- 
iment fiUl . we have a direct access to the carrier dynamics: Carrier 
distribution for different delay times after an optical excitation with 
a pump fluence of approximately 2.5 mJ/cm (marked blue in (b)). 
Gain is obtained during the first 300 fs for earner energies smaller 
than 250 meV. 



to the Dirac point resulting in a population inversion. The 
temporally subsequent decay of the population inversion on 
a hundred femtosecond timescale can be understood as an 
efficient interplay of Coulomb-induced Auger recombination 
1 121 UM (carrier loss) and intraband carrier-phonon scatter- 
ing (ener gy l oss). Signatures of these collinear scattering 
channels 1 1711 w ere recently also observed in photo-current 
experiments | Ijl and optical pump-probe measurements 

i nm. 



Applying the graphene Bloch-equations 123L 12411 . we mi- 
croscopically consider the light-carrier interaction as well as 
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carrier-carrier and carrier-phonon scattering in a consistent 
treatment 11251 12611 . A coupled set of differential equations 
describes the temporal evolution of (i) the occupation proba- 
bility p£ in the state k in the conduction and the valence band 
(A = c, v), (ii) the microscopic polarization reflecting the 
transition probability between both bands, and (iii) the phonon 
occupation with the momentum q for different modes j 

pi = 23 [nrnj + rL" A [1 - pi] - r k >£, 

pu = [tAoj u - i,ni x - 7k ] p k - tfll c [ Pk - pi] + Z4, (!) 

pern fi i ^3 1 pa63„ 
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The coupling to the exciting pulse is described by the Rabi 
frequency S! k c and its corresponding intraband contribution 
n£ x . For our analysis we generate a non-equilibrium carrier 
distribution vi a op tical excitation, similar to the experimen- 
tal realization 11311 : a laser pulse in $l£ A with a duration of 
10 fs and a photon energy of 1.5 cV is applied. A sketch of 
excitation and the linear dispersion is shown in Fig. QI a )- A 
delayed probe pulse (blue) is used to test the pump (red) in- 
duced population change. Many particle interactions (r, Uu, 
and 7k) are treated within second-order Born-Markov approx- 
imation, where we take into account a self-consistently deter- 
mined many-particle broadening of the strict energy conser- 
vation 12711 . We obtain a Boltzmann-like scattering equation, 
where the time- and momentum-dependent scattering rates 
r^™* consider carrier-carrier as well as carrier-phonon con- 
tributions. A similar equation is obtained for the phonon oc- 
cupation with the emission and absorption rates p^™/ abs The 
microscopic polarization pk is damped by the many-particle 
diagonal dephasing 7k and deformed by the off-diagonal con- 
tribution Wk 12311 . In exception of the linear transition fre- 
quency Awk, the experimentally estimated phonon decay rate 
Iph 12811 . and the Bose-Einstein distribution n J B as a thermal 
bath for the dynamic phonon occupation, all terms in Eqs. (Q~|i 
are explicitly time-dependent. A detailed description of these 
equations can be found in Ref. fl23ll . 

For the understanding of the formation of optical gain 
we analyze as an experimentally addressable observable the 
pump pulse-induced population inversion Apk(t) at the posi- 
tion of the probe field |ko | = ujl/2vf, where lol is the probe 
frequency and vp is the carrier velocity in graphene's linear 
band structure. From the numerical solution of the graphene 
Bloch equations (HJ we obtain Apk (i), which is approxi- 
mately proportional to the differential optical response, e.g. 
the reflection AR/R JUS. 

Figure |TJb) shows Apk (t) for different pump fluences up 
to 4mJ/cm 2 . In excellent agreement with the experiment 
1 1311 . we observe in the low excitation regime a linear scaling 
of the peak of Apu (t) with the pump fluence, and a clearly 
nonlinear saturation effect in the range of 0.5 — 1 pJ /cm 2 . 
At very high fluences, we find first signatures of Rabi oscilla- 
tions. The gain (population inversion) regime can be directly 
observed in the energy-resolved carrier distribution, shown 



in Fig. [TJc) for an exemplary pump fluence of 2.5 pj/cm 2 
(marked blue in Fig. \V[b)) at different times after the opti- 
cal excitation. At time t = 0, i.e. at the maximum of the 
excitation pulse, a peaked non-equilibrium population distri- 
bution is found. The low energy occupation wing stems from 
the initial thermal distribution at room temperature broadened 
by relaxation processes akeady acting during the pulse. In 
the course of time, after the pulse, the low energy population 
increases and we observe a population inversion for carrier en- 
ergies of up to 250 meV. An analysis of the distribution shows 
that as long as the population inversion is present no Fermi- 
Dirac distribution is established. A non-equilibrium between 
the carrier and phonon system is present. Depending on the 
pump fluence, it takes a few hundreds of fs until equilibrium 
is reached. 

Next we turn to the microscopic explanation of the forma- 
tion mechanism of the population inversion. A detailed in- 
spection shows, that the population inversion at lower energies 
is a direct consequence of energy loss of the optically excited 
carriers due to intraband carrier-phonon scattering, cf. Fig. 
|2a) (green, left). This process fills the states with low mo- 
menta having a low (linear) density of states. Hereby, the in- 
terplay of the reduced density of states (vanishing at |fc = 0) 
and the conserved carrier density results in a higher occupa- 
tion at lower energies in comparison to the optically excited 
distribution, where no population inversion is found. This 
process of energy loss within a band is indirectly accelerated 
by an ultrafast broadening of the non-equilibrium carrier dis- 
tribution via intraband Coulomb-scattering. A counteracting 
process is phonon-induced recombination between the bands, 
cf. Fig. |2J a ) (blue, middle). However, in graphene it is less 
efficient due to the reduced density of states close to the Dirac 
point, resulting in a phonon bottleneck 131 1. 

This overall picture of the population inversion mechanism 
is supported by Fig. [2] The population distribution 50 fs af- 
ter the excitation pulse under consideration of different re- 
laxation channels is shown in Fig. |2jb). For the purple line 
all scattering channels are considered and a population inver- 
sion is obtained for carrier energies smaller than 250 meV. 
Next, we neglect carrier-phonon scattering to demonstrate its 
crucial role for the build-up of a population inversion (blue 
dashed-dotted line in Fig. |2jb)). Without this mechanism of 
intraband energy loss, the population inversion is suppressed 
by efficient Coulomb-induced Auger recombination 132143411 . 
sketched in Fig. |3a) (red). Its redistribution of carriers to 
higher energies results in a single Fermi distribution for con- 
duction and valence band with a vanishing chemical potential 
1 3411 . cf. Fig. I2b). Taking into account all relaxation channels 
but the Auger processes (red dotted line) gives rise to an ener- 
getically broad gain regime with an unrealistically extensive 
persistence. This reflects the role of Auger-induced recombi- 
nation to the transience of the population inversion. 

One direct consequence of the carrier recombination 
(Auger and electron-phonon) in the process of transient popu- 
lation inversion, which is in full agreement with experimental 
measurements 11311 . is a decay of the population inversion on 
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FIG. 2. (a) Sketch of the most important scattering channels: en- 
ergy loss due to intraband phonon emission (green, left), phonon- 
induced recombination (blue, middle), and Coulomb-induced Auger 
recombination (red), (b) Carrier distribution 50 fs after the excitation 
pulse (2.5 /ij/cm ) considering the full dynamics (purple line), only 
Coulomb-induced scattering processes (dashed-dotted blue line), and 
all relaxation channels in exception of Auger-processes (dotted red 
line). The figure illustrates that the intraband carrier-phonon scat- 
tering is responsible for the occurrence of the population inversion, 
while Auger processes are counteracting, (c) Coulomb- (red) and 
phonon-induced (blue) recombination rates, reflecting their contri- 
bution to the decay of the gain regime, where the process of Auger 
recombination is the predominant mechanism. 



a femtosecond timescale. Such a decay has been measured 
in the pump-probe signal in Ref. 1 1311 . To study the decay 
of the population inversion in more detail, we analyze the 
phonon- and Coulomb-induced recombination rates, shown in 
Fig. |2jc). Our calculations reveal that in the strong excitation 
regime Coulomb-induced Auger processes (red line) are pre- 
dominant recombination channels, whereas interband carrier- 
phonon scattering (blue line) plays a minor. In particular, 
immediately after the excitation pulse the Coulomb-induced 
Auger recombination reaches values of up to 4 • 10 11 /fs cm 2 , 
which are large compared to the phonon-induced recombina- 
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FIG. 3. Maximal optical probe energy, where a population inversion 
can be obtained depending on the pump fluence. The values are ex- 
tracted from the carrier distribution 40 fs after the excitation pulse, 
corresponding to the temporal resolution of the experiment in Ref. 



tion of only 2 • 10 10 /fscm 2 and the overall carrier density 
of 4 • 10 13 /cm 2 . Consequently, the Auger recombination 
outweighs the phonon-induced intraband filling of the low- 
energy states resulting in a femtosecond decay of the popula- 
tion inversion. 

A key quest concerning applications of graphene in opto- 
electronics is the frequency regime where optical gain can be 
obtained and its dep endence on the excitation fluence. As re- 
ported in Ref. 11311 . for a certain excitation strength the car- 
rier occupation at a probe energy of 1.55 eV remains beneath 
the gain regime, whereas at lower probe energies a population 
inversion was detected. We use our access to the momentum- 
resolved carrier distribution to determine the maximal optical 
probe energy, at which population inversion is obtained de- 
pending on the pump fluence, shown in Fig. [3] The values 
are taken 40 fs after the excitation pulse, corresponding to the 
temporal resolution of the experiment lfl3ll . Up to a critical 
pump fluence of approximately 200/iJ/cm 2 we observe no 
population inversion. Above this value the population inver- 
sion appears abruptly at low probe energies. For stronger ex- 
citations the gain regime is broadened up to 0.75 eV. The ex- 
perimentally observed population inversion at a probe energy 
of 1.16 eV can be obtained for pump fluences > 5 mJ/cm 2 . 

Finally, we analyze the energy-resolved carrier distribu- 
tion during and after the optical excitation with an increas- 
ing pump fluence, cf. Fig. SJa-d). The shown occupations 
along the direction of the strongest li ght- matter coupling in 
the two-dimensional reciprocal space |23[I directly illustrate 
the action of the exciting laser pulse 13511 . The rise of the 
non-equilibrium occurs wave-like reflecting the optical pump 
frequency. In exception of panel (a), where a moderate excita- 
tion strength beneath the critical value was used, we observe 
a gain regime (red area) in a broad energy range depending 
on the applied pump fluence. Note that we also obtain popu- 
lation inversion at the excitation energy decaying on the time 
scale of the pump pulse duration. However, this should not be 
visible in experiments using the same temporal resolution for 
pump and probe pulse. Furthermore, the energy and time res- 
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FIG. 4. Time- and energy-resolved carrier distribution along the 
direction of the strongest light-matter coupling during excitations 
with different pump fluences: (a) moderate excitation regime, (b) 
gain regime, (c) Rabi-oscillation regime, and (d) second-harmonic 
regime. 



olution of the carrier distribution uncovers different aspects of 
the light-matter coupling in strong excitation regime. The ex- 
citation used in Fig. |4jc) gives rise to Rabi oscillations, since 
optical pumping overweights the diagonal dephasing 7k for 
the microscopic polarization , cf. Eqs. ((T). For the highest 
pump fiuence, shown in Fig. SJd), even a wing of the coherent 
two-photon absorption is visible. 

In conclusion, our calculations confirm the possibility of 
transient optical gain and population inversion in optically 
pumped graphene. We identify intraband energy loss via opti- 
cal phonons as the predominant mechanism for the occurrence 
of the population inversion. Its decay on a femtosecond time 
scale is mainly driven by the process of Auger recombination. 
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